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The designs of the first generation of cosmological 21-cm observatories are spht between single 
dipole experiments which integrate over a large patch of sky in order to find the global (spectral) 
signature of reionization, and interferometers with arcminute-scale angular resolution whose goal 
is to measure the 3D power spectrum of ionized regions during reionization. We examine whether 
intermediate scale instruments with complete Fourier (uv) coverage are capable of placing new 
constraints on reionization. We find that even without using a full power spectrum analysis, the 
global redshift of reionization, Zreion, can in principle be measured from the variance in the 21-cm 
signal among multiple beams as a function of frequency at a roughly 1 degree angular scale. At this 
scale, the beam-to-beam variance in the differential brightness temperature peaks when the average 
neutral fraction was ~ 50%, providing a convenient flag of Zrdon- We choose a low angular resolution 
of order 1° to exploit the physical size of the ionized regions and maximize the signal-to- noise ratio. 
Thermal noise, foregrounds, and instrumental effects should also be manageable at this angular 
scale, as long as the uv coverage is complete within the compact core required for low-resolution 
imaging. For example, we find that Zrcion can potentially be detected to within a redshift uncertainty 
of Azreion ^ 1 in ^ 500 hours of integration on the existing MWA prototype (with only 32 x 16 
dipoles), operating at an angular resolution of ~ 1° and a spectral resolution of 2.4 MHz. 



I. INTRODUCTION 

The 21-cni line from neutral hydrogen is a powerful 
probe of the high-rcdshift universe. Its observations 
could potentially constrain the properties of the first 
galaxies flWS], the cosmological initial conditions [MS], 
and possibly fundamental physics |71 |H| . There are cur- 
rently two standard approaches to the detection of the 
cosmological 21-cm signal: measuring the redshift evo- 
lution of its sky-averaged value (commonly labeled the 
"global 21-cm signal") 9 , and measuring the statistical 
properties of the full three-dimensional fluctuations in 
the 21-cm intensity The leading global experiment, 
EDGES, uses a single dipole to constrain the duration 
of reionization [TUHIS]- More advanced single-dipole ex- 
periments, such as a lunar satellite, are currently being 
planned fTT. Multi-dipole interferometers such as LO- 
FAR MWA [16 , GMRT [17], and PAPER [18], will 
search for the fluctuations caused by ionized regions dur- 
ing the epoch of reionization (EoR) and their 3D power 
spectrum [6lfT9ll20]. 

In both approaches, a central technical challenge is to 
remove the galactic synchrotron and other foregrounds 
from the cosmological 21-cm signal (CS) associated with 
the reionization of hydrogen |2r . The mean value of the 
foregrounds are on the order of hundreds of Kelvin even 
in quiet parts of the sky at ^ lOOMHz, while the magni- 
tude of the CS is on the order of tens of mK[Hl[Hllini- The 
foreground spectrum is smooth in frequency relative to 
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the CS, which makes foreground subtraction possible in 
principle [TT] [Ml [22 |2S] . Generically, the measured sig- 
nal can be fit by a slowly varying function in frequency 
(for instance, a polynomial in frequency or log frequency 
[TT[ [23], or non-parametric fits [5S]), with the residuals 
containing some combination of fitting error and the fiuc- 
tuating CS. If the fit is inadequate, residual foregrounds 
contaminate the CS. If the fit too closely matches the 
observed data, the CS will be removed together with the 
foreground. Several authors have already investigated 
the optimal fitting techniques that offer a good compro- 
mise for the tension between these opposing restrictions 

[13125]. 

So far, a single dipole antenna such as in EDGES con- 
strained the duration of reionization to be longer than 
Az > 0.07 |13 . The need to fit the unknown foreground 
and instrument response by a high-order polynomial in 
frequency suppresses the ability to detect a gradual reion- 
ization process as a function of redshift [51 [HI [23 ■ Only 
a sharp change in the mean brightness temperature as 
a function of frequency is detectable, making it difficult 
to constrain the most likely theoretical scenarios [TOl [26] . 
Thus, at present, we still do not have a 21-cm detection 
of the "redshift of reionization" - the point at which the 
universe was half ionized. 

Given the above experimental status, we would like 
to examine whether a modest interferometer (such as 
an early stage prototype of an interferometer like MWA 
32T) can do any better than a global experiment like 
EDGES at detecting the redshift of reionization. While 
an interferometer cannot measure the global signal, it 
could trace evolution in the statistical properties of the 
21-cm brightness fiuctuations with better calibration and 
foreground subtraction prospects than for a single dipole 



2 



experiment. In this paper, we calculate the strength of 
the 21-cm signal in a given fiducial model and give a 
plausibility argument that a modest (512 dipole) inter- 
ferometer with complete Fourier (uv) coverage should be 
adequate to measure the redshift of reionization in that 
scenario. Measuring Zicion would be a natural first tar- 
get for 21-cm surveys and an ideal way to demonstrate 
proof-of-concept for low-frequency interferometers at an 
early stage. 

We make several simplifications in our calculation. 
Most importantly, we assume full or nearly full uv cov- 
erage, which eliminates the influence of sidelobes. Fortu- 
nately, this assumption corresponds to the typical design 
goal of first generation low-frequency arrays being built 
for EoR science [17]. Under this assumption, our results 
do not depend fundamentally on the details of tile place- 
ment and rotation synthesis. As long as any particular 
array is designed to have a nearly complete uv coverage, 
our calculation should be applicable. It has also been 
shown in foreground and point source removal studies 
that full uv coverage greatly mitigates other instrumen- 
tal problems, which we also do not model explicitly. 

Since we are interested in full uu-coverage, all of the 
information collected by the interferometer is contained 
in the "dirty map" of the sky at the maximum resolu- 
tion of the instrument. Since a prototype instrument 
will not be able to probe down to the arcminute scales of 
reionization (where the variations are the strongest) and 
still maintain full uv coverage, most of the information 
of interest will be in the smallest scale available to us. 
The pixel-to-pixel variance of a "dirty map" is a simple 
way to encapsulate the information at the scale of inter- 
est [5S]. While this does throw out some information, 
for instance the location of the pixels and the longer- 
wavelength Fourier modes, it is a simple statistic to work 
with. If a detection can be shown to be feasible with 
just the variance of the dirty map, then it will surely be 
feasible with a more complete analysis. 

Specifically, our approach is to calculate the standard 
deviation of the residuals after foreground subtraction in 
each dirty map, as a function of frequency. What will 
remain in these residuals are the brightness temperature 
fluctuations (cTb) on top of a noisy spectrum and any 
fltting residuals. We show that in the regime where the 
spin temperature of the cosmic gas is much higher than 
the CMB temperature {Tg ^ Tcmb), a peak in axb is 
reached when the ionization fraction reaches 50%, as 
noted by others [57], providing a natural flag for Zicion- 
This holds when the pixels of the map are independent, 
which is roughly valid as long as their size is not less than 
the size of the ionized bubbles [53]. 

While ^rcion has been constrained by the measurement 
of the total optical depth for electron scattering with 
WMAP satellite, this is in fact an integrated constraint 
on the entire history of reionization [H |30]. The tech- 
nique considered here picks out a specific redshift with- 
out integrating along the line of sight. By combining the 
two measurements, it should be possible to constrain the 



duration of reionization as well. 

In order to estimate the prospects for this kind of mea- 
surement, we consider ID statistics representing the evo- 
lution of a pixel with frequency (a beam) . We investigate 
how these single-beam ID-statistics vary as a function of 
angular resolution (or "beam width") and consider in- 
termediate beam widths between the standard global ex- 
periments and the full MWA or LOFAR arrays. We also 
consider the dependence on frequency resolution. This is 
necessary because both the thermal noise and brightness 
temperature fiuctuations depend strongly on the angular 
scales involved. Signal fiuctuations increase with higher 
resolution, because the instrument averages over less of 
the sky. However, this dependence is not trivial because 
the bubbles have a finite size in real space and their loca- 
tions are correlated. Since the noise temperature varies 
according to the radiometer equation, thermal noise in- 
creases rapidly at higher resolutions at a fixed collecting 
area and integration time. To calculate the expected sig- 
nal, we adopt a semi-numerical realization of reionization 
and idealize 21-cm observations as pencil beams passing 
through the simulation box. 

The beam width dependence of arb has been previ- 
ously studied with the aid of a radiative transfer simula- 
tion in a cosmological box of 4/i~^ comoving Mpc (cMpc) 
on a side |31) . However, this size is an order of magni- 
tude too small to resolve the typical bubble sizes at the 
end of reionization, and is thus of limited applicability. 
Simulation boxes of at least 100 comoving Mpc are nec- 
essary to reliably capture the bubble size distribution at 
the end of reionization [32]. We use a publicly available 
semi- numerical code called 21cm FAST to overcome this 
hurdle [SSIIM]- 

The outline of this paper is as follows. In Sjlljwe out- 
line the methodology used as a function of beam width 
and frequency resolution. In §III| we consider the signal 
in a typical beam and calculate the effect of foreground 



subtraction on it. In {TV we explore the statistics of an 
ensemble of one-dimensional beams. In fjV] we explore 
how to optimize beam width and evaluate prospects for 
In SjVIl we briefly consider instrumen- 



measurmg Zj-oi 
tal issues and show that they should be surmountable 
in light of previous work. Finally, §VII| summarizes our 
main results and caveats. 



II. METHODOLOGY 

In order to study the statistical properties of the CS, 
we simulate a cosmological box using the publicly avail- 
able 21-cm FAST code [33 El]. This code generates a 
random realization of the matter density field and then 
uses an excursion-set approach to identify ionized re- 
gions. It has tunable cosmological parameters, ioniza- 
tion prescriptions, and several optional features such as 
a halo-finder and a spin temperature evolution tracker. 
We do not use these optional features and instead sim- 
ply use default cosmological and reionization prescription 
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settings, documented in Ref. (31]. For our analysis, we 
use two different realizations, one with merely a linear 
density evolution and the other with nonlinear evolution. 
Using both a linear and non-linear evolution enables us 
to resolve a large volume and achieve high resolutions for 
our simulated beams in the redshift range of 7 < z < 15. 

Comparisons to radiative transfer simulations have 
shown that this technique is reasonably accurate [331 El] ■ 
For our purposes, this method has the advantage of be- 
ing very computationally inexpensive while maintaining 
accuracy at the tens of percent level |35| . an advantage 
shared by this class of programs [3B^ . 

Our non-linear high resolution simulation involves a 
300 cMpc cosmological box with 600 pixels on a side, i.e. 
0.5 cMpc per pixel (corresponding to a frequency reso- 
lution of l^Vmin ~ 28 kHz, or an angular resolution of 
&min ~ 0.2' in this redshift range). Our larger, low reso- 
lution, linear realization contains 600 cMpc with a reso- 
lution of 200 pixels on a side, or 3 cMpc per pixel (cor- 
responding to Af,„i„ « 170kHz, or /S.9min ~ 1.1')- We 
then generate simulated sight-lines ("skewers") of vary- 
ing widths and line-of-sight depths which start at a ran- 
dom location and orientation. For simplicity, we adopt 
a square, top-hat cross-sectional (transverse to the line- 
of-sight) shape to our skewers and assume a number of 
beams for an instrument which has a 16 square degree 
field-of-view. 

We then assume randomized periodic boundary condi- 
tions on the box, such that if a skewer reaches the edge 
of the box, it re-enters the box at a new random orien- 
tation, face, and location. The redshift is advanced in 
discrete time steps of Az = 0.25 (corresponding to a fre- 
quency interval of 3.5 MHz at the redshifts of interest), 
and the sampling process resets after each time step. We 
sample each redshift in proportion to the number of real- 
space pixels which would be crossed during that redshift 
interval. Until the redshift advances a step, we linearly 
interpolate to calculate a frequency for that part of the 
beam. 

We assume that the beams have a fixed comoving 
width, each corresponding to a fixed angular scale at the 
high redshifts of interest (for z « 10, the variation in an- 
gular size for a constant comoving width between z — 1 
and z = 12 is < 25%). We adopt ACDM parameters of 
erg = 0.817, h = 0.7, an = 0.28, ^ ^ 0.046, 17a = 0.72, 
consistent with the 7- year WMAP data [30] . Throughout 
the paper, all length scales are in comoving units unless 
otherwise specified. 



III. SINGLE-BEAM PROPERTIES 

Each observing beam from an interferometer samples a 
new region of space. Therefore, one should expect cosmo- 
logical fluctuations in the differential brightness temper- 
ature within each beam in frequency space. The 21-cm 
line is observed in either emission or absorption, depend- 
ing on whether the local spin temperature, Tg, is above 



or below the cosmic microwave background (CMB) tem- 
perature, = 2.73K X (1 -I- z). The 21-cm fluctuations 
corresponding to peculiar velocities, density fluctuations, 
and ionization fraction fluctuations are described by the 
familiar equation, ^ 

5T,« 9x^,(1 + '5)(l + z)V2 

( H{z)/{l + z) \ ^ ^ 

n dv,/dr, ) 

where STjj is the differential brightness temperature, 
xhi = ^—Xi is the neutral hydrogen fraction (by volume), 
5 is the overdensity, and H{z) is the Hubble parameter 
at redshift z which is equal to the line-of-sight velocity 
gradient (1 -f z)dv\\/dr\\ in the absence of peculiar veloc- 
ities [S]. In the linear regime applicable on >Mpc scales, 
density fluctuations can be ignored relative to the order 
unity fluctuations in the neutral fraction xhi- If peculiar 
velocities are also ignored and Tg » T^, then 

5Tb^QxHi{l^ zf'^ mK. (2) 

The brightness temperature fluctuations over scales of > 
Mpc can be modeled as overlapping, growing spherical 
"bubbles" of ionized hydrogen which form around the 
flrst ionizing sources [8l fTOl l22l [37 l l38] . 

The 21-cm signal from reionization gives each beam 
a complex structure, because each ionized bubble can 
occupy part or all of the beam both along the line of 
sight and in the plane of the sky. In addition, overlapping 
bubbles might correlate JTf, among adjacent beams which 
belong to the same bubble [20] • This correlation can be 
characterized by the form |38) : 

(x,(r)x,(r')) =xj + {x, - xj)f(r/R,), (3) 

where r = |r — r'\ and f{r/Rc) parametrizes the effect of 
the flnite characteristic size Rc of the bubbles on the cor- 
relation. The characteristic scale of the H II bubbles has 
been shown to depend primarily on the mean ionization 
fraction and, to a lesser degree, on other parameters |22J. 
Equation [3] implies that if r 3> Rc, f{^/Rc) approaches 
zero and the variance of STj, is maximized when the ion- 
ized fraction is 50% (ignoring density fluctuations). We 
will use this feature as our marker for the reionization 
redshift, Zreion, and verify how valid this is at the rele- 
vant scales using 21-cm FAST. 

Figure [l] shows typical skewers (beams) with four dif- 
ferent widths. Two effects are evident from these ex- 
amples. First, the underlying bubble structure causes 
the CS to vary with frequency much more rapidly than 
its mean behavior, even though reionization itself is pro- 
gressing slowly. Secondly, it is visually obvious that while 
the mean behavior of these signals is the same, the mag- 
nitude of the variations from the mean depends strongly 
on beam width, which we consider more carefully in §IV| 

In order to recover these skewers, it is necessary to sub- 
tract out the foregrounds using some kind of flt or model 
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of them. While others have done this in much greater 
detail, we briefly verify that this is possible in princi- 
ple in our analysis - i.e. even a high-order polynomial 
will not remove all the power from the CS fluctuations 
in frequency space. In order to see what part of the sig- 
nal would remain after a foreground subtraction, one can 
add an approximation of the noise, and then fit the com- 
bined signal and subtract out the fit. We approximate 
the signal in quiet portions of the high-latitude radio sky 
as a power law in frequency [8]: 
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After adding this mean behavior to our individual beams, 
we subtract out a 12th order polynomial in frequency to 
fit the signal, and plot the residuals in Figure [2j 
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FIG. 1. Typical skewers with beam widths 6^ of 1' (top 
left), 9' (top right), 1.2° (bottom left), and 2.4° (bottom 
right), with A;/ ~ 30kHz. Increased beam width reduces 
the average level at which the signal varies by averaging over 
a bigger area of sky at each frequency interval. There is no 
thermal or instrumental noise in these examples. As others 
have found, foreground subtraction on these signals with high 
order polynomials does not remove much of the power from 
the line-of-sight fluctuations. 

It is evident from Figure [2] that the fitting residuals 
do not in principle dominate the signal if all we are con- 
cerned about is the magnitude of the sky temperature 
varying with frequency. We did not make a fully ade- 
quate foreground model here, but we find fitting residuals 
to be at about « 0.75 mK for the most relevant angular 
resolution of 1.2° with a 3rd order polynomial, or k, 0.4 
mK with a 3rd and 12th order polynomial, with fore- 
ground subtraction performed before binning of adjacent 
frequency bins. 
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FIG. 2. Signal after subtracting out a 12 order polynomial 
fit to the spectral data as a function of beam width. In each 
panel, the dashed top curve shows the original signal with- 
out any foregrounds, while the bottom solid curve shows the 
residuals after foreground subtraction. The beam widths at 
z = 9 are 0.3' (top left), 3.4' (top right), 10' (bottom left), and 
0.6° (bottom right). Note the large difference in the residual 
signal amplitude between the different panels. 



Of more concern is whether foreground subtraction is 
still effective after including point sources and other fre- 
quency dependent distortions. This question is studied 
much more carefully by other authors |391 140j in the con- 
text of a full first generation array. With 3rd order poly- 
nomials. Bowman et al. find a best case scenario of 1 mK 
residuals after foreground subtraction. The magnitude 
of these residuals sets a systematic bound on the min- 
imum signal strength that will be distinguishable from 
foreground subtraction errors (as discussed in SV]). We 
take the approach of trying to keep the magnitude of the 
signal above 1 mK while noting that increasing the order 
of the polynomial could possibly mitigate this problem 
(with caveats about losing more of the long-wavelength 
signal) . 

It remains to be seen whether a prototype array will 
be able to overcome these issues, although we are hopeful 
that it will (see ^ VI). At any rate, our results demonstrate 



that very little power from the CS (from the scales of in- 
terest to us) is removed within a single beam, even by 
a high-order polynomial at degree-scale resolution. Fig- 
ure [2] provides a simple way to understand why global 
experiments have a difficult time detecting reionization 
in a gradual reionization case: as the resolution scale in- 
creases, the fiuctuations get smaller, making it more dif- 
ficult to separate the post-subtraction fluctuations from 
noise alone. 
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Thus, our post-subtraction signal is composed of these 
fluctuations plus thermal noise in the antenna which also 
generates deviations. To demonstrate that we can mea- 
sure these fluctuations accurately, we need to quantify 
the thermal noise in the antenna due to the mean sky 
temperature and the magnitude of deviations as a func- 
tion of beam width and frequency resolution. We address 
the latter issue flrst. 



IV. BEAM-TO-BEAM VARIANCE 

Next, we consider an ensemble of beams that are suf- 
flciently narrow to probe fluctuations caused by inhomo- 
geneous reionization. We can then address the "beam-to- 
beam" variance (u^i,), namely the variance in the signal 
at a given frequency bin over an ensemble of beams. This 
is a natural quantity to consider for an interferometer 
which observes a number of beams, A^bcams ~ ^^/(^w: 
where Afl is the solid angle of the field of view and 6^ 
is the width of each square beam. The beam-to-beam 
variance or standard deviation is useful because inter- 
ferometers only measure the deviations from the mean 
signal. 

Since the number of regions averaged over is a function 
of the resolution, it is important to consider the beam 
width and frequency dependence of the beam-to-beam 
variance. If there were no correlations among nearby 
regions (i.e., bubbles were much smaller than the re- 
gions under consideration, and regions were not corre- 
lated), the sample standard deviation would fall off as 
1 / \/ ^pixels where iVpixeis is the number of pixels. How- 
ever, inside of large ionized regions, there is near-perfect 
correlation among nearby pixels. Moreover, because the 
ionized regions themselves are correlated, the maximum 
variance may occur in an intermediate regime, where just 
a few bubbles are being averaged over. The largest ion- 
ized bubbles obtain a size of ^ 100 cMpc at the end of 
reionization [41]. This corresponds to ^ 0.6 degrees at 
z=9. Thus, the exact dependence of axb on resolution is 
both non-trivial and important to this measurement. It 
is most convenient to use 21-cm FAST to estimate this 
standard deviation of (5Tf, as a function of both frequency 
and angular scale. An illustration of this dependence is 
depicted in Figure [Sj 

As expected, the beam-to-beam standard deviation de- 
creases more slowly than the 1 / iVpixeis statistics would 
predict. For instance, in the bottom right panel of Fig- 
ure |3j it takes a factor of ~ 8 more pixels to reduce the 
signal by a factor of ~ 0.5. Each pixel is being randomly 
chosen, but clusters of pixels chosen together will be cor- 
related because of the finite size of the ionized bubbles. 
As shown in the next section, this significantly boosts 
our signal-to-noise ratio, because the experiment requires 
seeing the peak of this signal against fluctuations which 
do obey 1 / -^/A^pixeis statistics. 

The location of the peaks on the right hand side of 
Figure [3] does seem to be moderately scale-dependent. 
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FIG. 3. Top left: Evolution of the mean ST^ in a 600 cMpc 
box. Top right: The standard deviation of STb as a func- 
tion of redshift for 6 beam widths (top to bottom: 1', 2', 
4', 0.6°, 1°, 2°). The frequency resolution was ~ 150 kHz in 
all cases. Note that for larger beam widths, the maximum 
variance occurs at a redshift of about 10 when the ionization 
fraction is 0.5. Bottom left: Evolution of the mean neutral 
fraction in a 600 cMpc box. Bottom right: The standard 
deviation of 5Tb as a function of redshift for 4 different fre- 
quency resolutions A;/ (top to bottom: 150 kHz, 300 kHz, 600 
kHz, 1.2 MHz). The beam width was 1' in all cases. Again, 
the peak occurs at z ~ 10 for large Ai/. 

The large component of the variance which originates 
from uncorrelated bubbles peaks at a neutral fraction 
of 50%. However, on any scale, there is an additional 
component originating from the second term of equa- 
tion ([3]) , which determines the redshift at which the vari- 
ance peaks. The inferred redshift of reionization (defined 
by the neutral fraction being 50%) has therefore some 
model-dependent uncertainty. For any specific model of 
reionization, this uncertainty can be removed by a com- 
parison to the model's version of Figure [3j 

V. MEASURING THE SIGNAL AGAINST 
THERMAL NOISE 

As mentioned above, our approach is to measure the 
standard deviation of each dirty map as a function of 
frequency. The standard deviation of the ensemble of 
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observed beams includes a few different components: the 
thermal noise, the standard deviation of 6T1, which peaks 
at some redshift, and any fitting errors, which we ig- 
nore. The variances of noise and signal at any redshift 
add directly by the Bienayme formula (since the noise is 
assumed to be uncorrelated with the signal): 

'AT^eas(f^) + systematic (5) 

The standard deviations of STb {(JTb) were calculated in 
[TV and are expected to peak near Xi = 0.5, provid- 
ing a flag of Zrcion- This peak sits on top of the noise 
spectrum, which should follow Eqs. ^ and ([6]). In an 
interferometer imaging experiment, the noise is obtained 
from the radiometer equation. Using typical values for 
the experimental parameters of 21-cm observations, the 
noise temperature is given by, [S]: 
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(6) 



where Atot is the effective collecting area of the inter- 
ferometric array, 9^ is the angular resolution, Av is the 
bandwidth, and ti„t is the integration time. At a fixed 
collecting area and integration time, the noise tempera- 
ture is very sensitive to the angular resolution and Av. In 
order to extract the standard deviation excess caused by 
the 21-cm fiuctuations from (Tdirty(i^): one can subtract 
off a model or power law fit to the {AT^){i^) spectrum. 
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We will assume that this fit will only leave residuals due 
to the finite sampling errors of measuring a variance. 
These sampling errors are the "effective noise" AT^s, 
which might cause us to make a mistake in determing the 
true variance at any frequency bin. For a Gaussian prob- 
ability distribution in the iVboams ^ 1 limit, the average 
deviation between the measured noise standard deviation 
and the true variance of the noise is given by a standard 
statistics result l42l: 
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Similarly, the sample variance of the signal would be: 
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(8) 
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Therefore, we can define the effective noise as 



(9) 



AT, 



eff 



1 



2iVbc 



-{{AT^Y + al,) (10) 



Under these assumptions, we define the effective signal- 
to-noise ratio as 



S/N, 
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^{{AT^f 
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(11) 



Our procedure is illustrated in Figure [4] We infer 
the redshift of reionization from the highest signal-to- 
effective-noise point on the curve, which, as mentioned 
before, is a somewhat model-dependent procedure that 
could be improved given some model for reionization. 
However, it assumes only the existence of a peak in the 
Ts 3> Tfc regime. In practice, systematics are likely to be 
the actual limiting factor in the choice of the resolution. 
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FIG. 4. An illustration of the different components in the 



measurement of the reionization redshift 



for the mod- 



est array parameters: Atot ~ 500 m^ (roughly 32x 16 dipoles), 
6^ = 1.2°, Ai/ = 2.4MHz, and Unt = 500 hours. The beam- 
to-beam rms variation (solid black top curve) is composed of 
both the variance due to thermal noise (dotted blue curve) 
and the signal caused by the bubbles (black short-dashed 
curve). If the expected variance from noise is subtracted off 
from the measured variance, all that remains are the signal 
and sampling errors as fitting residuals of the noise profile 
(red long-dashed curve). If there were no signal at all, the 
expected profile final result would only be due to sampling 
errors (bottom green dot-dashed curve). The signal-to- noise 
(purple dot-dashed curve) is calculated as defined in the text 
at each frequency bin. These noise estimates assume a "quiet" 
high latitude portion of the sky described by equation |4] 

We note that a higher angular resolution (small beam 
width 6w) increases the number of beams and the ab- 
solute magnitude of the signal. However, the signal-to- 
noise ratio increases at larger values of 9^ and Av up to a 
point, because the effective noise goes down more rapidly 
than the signal standard deviation up to the character- 
istic size of an individual ionized bubble. Furthermore, 
the scale must be big enough for the beams to be at least 
roughly independent, or the peak will not occur when 
the neutral fraction has reached 50%. Thus, to measure 
^roion, it is ideal to probe degree scales, at which signal- 
to-effective-noise-ratio is high, the beams are roughly in- 
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dependent, and the signal is comparable to the absolute 
magnitude of the noise and larger than systematic effects. 
In practice, we find this to be between ~ 0.6° and ^ 1.2° 
at Aiy = 2.4 MHz where the beams are roughly indepen- 
dent, (Tt6 is between 3 mK and 1 mK, and the signal 
to noise at the peak is between 11 and 13. Systematics 
will most likely limit us because foreground subtraction 
with low order polynomials most likely cannot be accom- 
plished at the better than 1 mK level [33]. If a stronger 
absolute magnitude of the signal is required, smaller an- 
gular resolutions can be used and compensated for by 
longer integration times. 

To estimate the error in Zj-cion, we randomly draw 
Gaussian noise, add it to the sample realizations shown 
in Figure [4j and calculate the standard error in the loca- 
tion of the highest peak after 100 Monte Carlo trials of 
this procedure. Our result is that for tint = 500 hours, 
Atot = 500 m^ (corresponding to MWA with 32 tiles of 
16 dipoles each), 6^ = 1.2°, and Ai' — 2.4 MHz, one can 
achieve a system noise temperature of AT^ — 0.94 mK 
and AT^ = 0.12mK at Zrcion- With these parameters, 
one can keep the CS above 1 mK and make a 4cr de- 
tection of Zrcion with a standard error of A^reion = 0.25. 
Depending on realistic integration time constraints and 
systematic issues caused by antenna layout, optimal pa- 
rameters may vary for actual experiments. The above 
numbers indicate that such a measurement is possible 
with current technology. However, these numbers ignore 
systematics from fitting errors, any residual discrepancy 
between the peak and the redshift of reionization, and 
RFI. 

A plot of the signal-to-noise ratio and absolute am- 
plitude of the signal, max(crTb) , versus increasing beam 
width and Av is shown in Figure [5] In Figure [6j we show 
how strong the constraint is on Azroion as a function of 
integration time for the observing parameters we chose. 



VI. PRACTICAL CONSIDERATIONS 

So far, our estimates have been restricted to mostly 
theoretical considerations about signal and thermal 
noise. In practice, however, detailed foreground and in- 
strumental considerations are critical to determining the 
plausibility of such a measurement of ^^rcion- 

The chief analytical technique which will reduce these 
effects is the standard signal fitting and subtraction tech- 
nique. Each pixel can be independently fitted with a 
suitable function (polynomial or otherwise) as a function 
of frequency, and contributions from unresolved point 
sources, galactic synchrotron emission, and frequency- 
dependent instrumental response should be reduced, 
leaving only the fast- varying CS [TTl [201 . 

However, it is also conceivable that a frequency- 
dependent point spread function (PSF) or field-of- 
view (FOV) could mix power from nearby pixels in a 
frequency-dependent way, and create artificial peaks in 
a measurement of axb versus frequency. This is of par- 



Ai/ (MHz) 
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T , ^ I ^ ^ ^ I ^ ^ , I ^ ^ ^ I ^ ^ ^ I ^ , ^ |— 
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I ' ' ' ' ' ' ' ' ' ' 1 

0.6 0.8 1 1.2 

0^ (degrees) 

FIG. 5. Estimated signal-to-noise ratio (dotted line) and ab- 
solute magnitude of signal (solid line) vs. beam width and 
band width. The angular width is increased along with fre- 
quency width to keep the ratio of the real space size of the 
pixels the same in different directions. The expected signal- 
to-noise ratio that can be achieved is set by the systematics of 
the experiment (left axis), as the absolute magnitude of axb 
needs to stay above the instrument systematic uncertainties 
(involving, e.g., residuals from foreground subtraction (right 
axis). In accordance with the radiometer equation, we assume 
that these limits will allow a signal of greater than 1 mK to be 
detectable. We do not show values of beam width for which 
we do not expect the beams to be independent, because this 
may introduce a non- negligible sampling error in art- 



ticular concern if the uv plane is sparsely filled, which 
would make the PSF deviate strongly from a S function. 
Also, if the uv plane is unfilled, pixels which may not 
be physically correlated due to the bubbles may become 
correlated due to under sampling. It is impossible to 
have more independent pixels than there are independent 
baselines, whatever the nominal angular resolution. 

Fortunately, we can rely on past work to show that 
this is unlikely to be an issue even for smaller prototype 
experiments. Bowman et al. found that, for a full MWA, 
a third order polynomial does an excellent job of sub- 
tracting foreground contamination for baselines where 
the uv coverage was complete [3^. Bowman et al. also 
found that for these larger angular scales, with uniform 
weighting of visibilities, the PSF was efi^ectively a clean 6 
function with no sidelobe confusion. Similar results were 
found by Liu et al. [iP] . 

In particular, using the formalism from Bowman et al. 
|39j , we can estimate the variance due to sidelobes in the 
synthesized array beam 

al^aUl + B'„ns^p/nB) (12) 
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FIG. 6. The standard error in the global redshift of reioniza- 
tion, Azrcion, as a function of total noise temperature at an 
observing frequency of 158 MHz, with 9^ = 1.2° and Au — 2.4 
MHz for Atot ~ 500m^ . The integration time in hours is given 
on the top axis. A^rcion is calculated by taking the standard 
deviation of fifty Monte Carlo trials, following the procedure 
illustrated in Figure |4] These errors ignore the systematic 
considerations discussed in the text. 

where Brms is the rms value of the beam response relative 
to the peak, fip corresponds to the FOV and ^ Gd 
is the solid angle of a resolution element. For MWA 32T, 
the value of Brms is ^ 1%, and Bl^^Vlp /VIb ^ 10% 
[43] , This implies the sidelobe confusion will be 10% of 
the sky variance, so sidelobe rms will be 5% of the sky 
rms, which is about lOK. In other words, the sidelobe 
rms will be '-^ 500mK, which is a reasonable worst-case 
scenario, and in practice one should be able to do better 
once foreground subtraction is taken into account. 

Therefore, one could hope an instrument with good uv 
coverage at large angular scales should be able to defeat 
these instrumental concerns. This is not a merely hypo- 
thetical scenario, as the uv plane should be reasonably 
well filled (about 60%) for a reference configuration of 
MWA 32T gg, as illustrated in Figure [t) Thus, it seems 
plausible that MWA 32T itself may have a good chance 
of detecting this signal, and an instrument with even bet- 
ter UV coverage would be more ideal. Perhaps LOFAR 
uv coverage |40| could be tuned to a similar performance 
level if long baselines were masked from the analysis to 
keep the uv plane filled. 

Since this approach has not yet been tested in practice 
there could be unknown instrumental issues, such as cali- 
bration issues for 32T, a non-flat bandpass (which would 
require a higher order polynomial), or issues with the 
beam response pattern i^B . There is certainly room for 
more study of these issues from an instrumental point- 



of-view. But it seems reasonable to hope that dense uv 
coverage of a compact core will ensure that frequency- 
dependent instrumental response will not hinder polyno- 
mial foreground subtraction. 



VII. CONCLUSIONS 

We have demonstrated in Figures |4] [5] and [6] that a 
modest interferometer with an optimal beam width and 
frequency resolution can statistically detect the global 
redshift of reionization, Zrcion- The standard deviation 
of dirty maps as a function of frequency at 1° scales 
shows that detection prospects are plausible and that 
a conventional power-spectrum analysis could probably 
do even better. In a single-dipole experiment, current 
foreground removal strategies for global 21-cm signal re- 
move much of the global signal itself unless reionization 
occurred over an unusually narrow redshift interval. Al- 
though the global signal has many interesting properties 
other than just the redshift of reionization, it can only be 
used to rule out fast reionization signal with the EDGES 
experiment (given current calibration). First-generation 
interferometers capable of measuring the 21-cm power 
spectrum in all its detail may still be several years away, 
but a detection may be possible now. 

We find that prototypes of interferometers designed to 
measure the power spectrum should be able to measure 
the global redshift of reionization without a large invest- 
ment of time or funds, based on the technique outlined 
in this paper. We quantified the observational prospects 
in a particular reionization scenario and showed that in- 
termediate beam widths of ^ 1.2°, integrated over hun- 
dreds of hours at MHz resolutions and ~ 500 m^ of ef- 
fective area should allow for a statistical detection of the 
global redshift of reionization if instrumental and system- 
atic challenges are overcome. 

It is difficult to predict how instrumental challenges, 
such as removal of point sources, uneven instrumental 
response, RFI, and other issues would play out for the 
measurement described here. But work done by others 
suggests that dense uv coverage with a compact core will 
mitigate many potential issues, and that foreground sub- 
traction should work at least as effectively as it does for 
an EDGES-like global experiment. 

Although the power spectrum of 21-cm fluctuations at 
all scales provides a more detailed probe of reionization, 
the ID flag of Zrcion considered here represents a "zeroth 
order" way to see that even prototype scale instruments 
can constrain reionization. In principle, one might be 
able to use the width of the axb peak to constrain 
the duration of reionization. This would be a simple 
way to demonstrate the feasibility of the interferomet- 
ric approach to the emerging frontier of 21-cm cosmology. 
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FIG. 7. Simulated uv coverage after rotation synthesis for 
MWA 32T. The upper panel shows the configuration (in x-y 
coordinates) of MWA 32T tiles, and the middle panel shows 
a snapshot of the corresponding baselines. The lower panel 
illustrates the tracks of the dipoles (due to Earth rotation) in 
the Fourier u-v coordinates, and illustrate that ~ 60% of the 
uv plane is filled after rotation synthesis with this prototype. 
The details of the attenna layout are not used in the present 
calculation. Plots are courtesy of J. Bowman |43| . 
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